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ABSTRACT 

This  paper  demonstrates  a  new  gravest  empirical  mode  (GEM)  technique  that  constructs  multi-index 
lookup  tables  of  temperature  ( T )  and  specific  volume  anomalies  (S)  using  historical  hydrocasts  as  a  function 
of  three  indices:  round-trip  travel  time  (r)  from  sea  floor  to  the  surface,  sea  surface  temperature,  and 
pressure.  Moreover,  the  historical  hydrocasts  are  separated  into  non-mixed-layer  (NML)  and  mixed-layer 
(ML)  groups,  and  a  single  GEM  field  is  constructed  for  each  group.  This  is  called  the  MI-GEM  technique. 
The  appropriate  dates  for  MI-GEM  fields  are  determined  by  the  monthly  distribution  of  the  number  of 
NML  and  ML  profiles  in  the  historical  hydrocasts,  which  are  also  well  correlated  with  the  strength  of  the 
winds  during  the  2  yr  of  observations.  The  T  and  S  profiles  that  are  determined  by  this  MI-GEM  technique 
capture  92%  and  88%  of  the  T  and  8  variances  in  the  depth  range  of  0-200  db.  These  values  reduce  by  about 
one-third  of  the  unexplained  error  variance  of  the  residual  GEM,  which  was  recently  developed  and  applied 
to  the  optimal  interpolated  t  data  in  the  southwestern  Japan/East  Sea  (JES)  by  Mitchell  et  al.  Comparisons 
with  the  in  situ  CTD  casts  demonstrate  that  the  MI-GEM  technique  almost  always  produces  improved  full 
water  column  profiles  of  T  and  S.  Whereas  the  residual  GEM  estimates  had  exhibited  qualitatively  erro¬ 
neous  features  like  T  inversions  in  the  near-surface  layer  and  too  thin  or  thick  intermediate  water  layers  in 
some  regions,  the  MI-GEM  estimates  avoid  those  problems,  which  were  inherent  to  the  residual  GEM 
technique  in  the  southwestern  JES. 


1.  Introduction 

The  gravest  empirical  mode  (GEM)  technique  is  a 
method  that  is  used  to  estimate  vertical  profiles  of  tem¬ 
perature  (T),  salinity  (5),  and  specific  volume  anomaly 
(8)  from  a  vertically  integrated  proxy,  such  as  round- 
trip  acoustic  travel  time  (t)  from  the  bottom  to  the 
surface,  geopotential  height,  or  heat  content  (Meinen 
and  Watts  2000;  Watts  et  al.  2001).  The  standard  GEM 
technique  produces  two-index  lookup  tables  (GEM 
fields)  of  T,  S,  and  8  from  historical  hydrocast  data  as  a 
function  of  t  (or  another  index  or  proxy)  and  pressure. 
This  standard  GEM  technique  has  been  applied  suc¬ 
cessfully  to  the  North  Atlantic  (Meinen  and  Watts 
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2000)  and  Antarctic  Circumpolar  (Sun  and  Watts  2001) 
Currents.  An  improved  standard  GEM  technique, 
which  removes  the  seasonal  signals  in  the  upper  layer 
prior  to  the  GEM  field  computation  to  reduce  the  as¬ 
sociated  errors,  has  been  developed  and  applied  to  the 
subantarctic  front  (Watts  et  al.  2001)  and  the  Kuroshio 
south  of  Japan  (Book  et  al.  2002).  However,  Mitchell  et 
al.  (2004)  recently  demonstrated  that  the  southwest¬ 
ern  Japan/East  Sea  (JES),  called  the  Ulleung  Basin 
(UB),  is  not  suitable  to  this  improved  standard  GEM 
technique  because  variations  in  the  shallow  permanent 
thermocline  occur  at  the  same  depths  as  the  seasonal 
signal,  which  varies  in  vertical  structure,  amplitude,  and 
phase  from  one  location  to  another,  even  within  this 
limited  region.  Therefore,  they  modified  the  method  to 
estimate  anomalies  relative  to  a  climatology,  referred 
to  as  the  residual  GEM  technique,  as  follows. 

The  residual  GEM  technique  builds  upon  a  climatol¬ 
ogy,  such  as  the  U.S.  Navy’s  Modular  Ocean  Data  As- 
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similation  System  (MOD AS).  Mitchell  et  al.  (2005) 
used  this  to  represent  the  spatiotemporal  annual  signal, 
and  used  historical  hydrocasts  to  estimate  the  residual 
(i.e.,  anomaly)  T  and  8  variability  from  those,  indexed 
by  residual  acoustic  travel  time  as  a  proxy  variable.  The 
residual  GEM  fields  were  additionally  grouped  under 
three  classes  of  residual  sea  surface  temperature  (SST), 
which  significantly  reduced  the  errors  near  the  surface. 
Mitchell  et  al.  (2005)  produced  accurate  daily  T maps  at 
a  100-m  depth  in  the  UB  for  2  yr,  using  the  residual 
GEM  technique  to  interpret  t  data  from  pressure  sen¬ 
sor-equipped  inverted  echo  sounders  (PIESs). 

In  the  UB,  dynamics  of  the  intermediate  water  is  of 
great  interest.  In  particular,  a  water  mass  with  a  vertical 
S  minimum  that  is  accompanied  by  a  dissolved  oxygen 
maximum  may  be  found  in  a  depth  range  of  100-400  m. 
This  water  mass  has  the  temperature  and  potential  den¬ 
sity  range  from  5°  to  1°C  and  from  26.9  to  27.3  ae,  and 
is  hypothesized  to  originate  from  subducted  water 
along  the  northern  region  of  the  subpolar  front.  Kim 
and  Chung  (1984)  found  it  in  the  UB  and  named  it  the 
East  Sea  Intermediate  Water  (ESIW).  Although  there 
are  several  observational  studies  showing  the  distribu¬ 
tion  and  circulation  of  the  ESIW  in  the  UB  (e.g.,  Kim 
et  al.  1991;  Shin  et  al.  1998),  those  hydrographic  studies 
produced  intermittent  snapshots  of  the  ESIW  circula¬ 
tion  that  were  insufficient  to  understand  its  source  and 
circulation.  In  contrast  to  the  above  studies,  the  dy¬ 
namic  height  estimates  from  our  PIES  array  can  pro¬ 
vide  daily  maps  of  the  circulation  at  intermediate 
depths.  However,  some  interpretations  that  used  the 
residual  GEM  estimates  were  found  to  be  erroneous 
regarding  intermediate  water  temperature  structures, 
by  estimating  an  unrealistically  thin  or  thick  ESIW 
layer.  Comparisons  of  the  residual  GEM  results  against 
concurrent  in  situ  hydrocasts  revealed  that  the  residual 
GEM  sometimes  produced  erroneous  T  and  8  inver¬ 
sions  near  the  surface  layer.  In  addition,  the  base  of  the 
thermocline  (pycnocline)  at  1°C  (27.3  cr0)  was  qualita¬ 
tively  too  shallow  or  too  deep  in  some  regions.  Al¬ 
though  the  residual  GEM  had  improved  upon  the  pre¬ 
vious  GEM  estimates  of  the  thermal  structure  in  the 
UB,  and,  indeed,  worked  quite  well  for  mapping  the 
temperature  at  a  100-m  depth,  it  unavoidably  produced 
the  inversions  and  other  errors  in  near-surface  waters 
and  was  complicated  to  apply. 

The  problems  with  the  residual  GEM  estimates  mo¬ 
tivated  us  to  develop  a  further-improved  GEM  tech¬ 
nique  that  is  applicable  to  the  UB.  In  this  paper,  we 
demonstrate  a  multi-index  GEM  technique,  which  uses 
satellite-mapped  SST  as  an  additional  parameter  to 
make  full  multi-index  lookup  tables  for  T  and  8  versus 
pressure,  and  also  distinguishes  between  non-mixed- 


layer  (NML)  and  mixed-layer  (ML)  profiles,  as  de¬ 
scribed  later.  Herein  we  call  this  the  MI-GEM  tech¬ 
nique.  We  apply  this  new  GEM  technique  to  the  t  mea¬ 
surements  from  an  array  of  PIES  in  the  UB. 

2.  Data  and  methods 

We  use  t  records  measured  by  an  array  of  PIES  in 
the  UB,  a  set  of  2248  historical  hydrocasts,  441  in  situ 
CTD  casts  from  the  Korea  Ocean  Research  and  Devel¬ 
opment  Institute  (KORDI)  and  the  Korea  Oceano¬ 
graphic  Data  Center  (KODC),  and  SST  from  the 
Tropical  Rainfall  Measuring  Mission  ( TRMM )  satellite 
observations.  These  data  are  the  same  as  those  used  in 
Mitchell  et  al.  (2004).  In  addition,  time  series  of  Navy 
Operational  Global  Atmospheric  Prediction  System 
(NOGAPS)  wind  stress,  on  a  degree  grid  with  a  time 
interval  of  12  h,  were  used  to  distinguish  time  intervals 
when  a  deepened  ML  should  occur.  A  brief  description 
of  t  and  SST  datasets  is  presented  here.  For  more  de¬ 
tails,  see  Mitchell  et  al.  (2004). 

An  array  of  25  PIESs,  designed  to  have  55-60-km 
spacing  covering  the  UB,  was  deployed  during  1999- 
2001  (Fig.  1).  The  PIES  measured  t  each  hour  with  a 
0.05-ms  precision.  Two  PIESs  at  P32  and  P41  were  lost, 
probably  because  of  deep-crab-fishing  activity.  Five- 
day  low-pass  filtering  was  applied  to  the  t  records  to 
focus  on  geostrophic  signals.  Then,  the  t  records  were 
subsampled  at  1200  UTC  every  day  for  736  days  from 
16  June  1999  to  20  June  2001. 

Application  of  optimal  interpolation  (OI)  mapping 
to  the  subsampled  r  requires  the  selection  of  a  stan¬ 
dardized  reference  level,  because  each  PIES  measured 
t  at  its  individual  bottom  depth.  Mitchell  et  al.  (2004) 
chose  a  reference  level  of  500  db  and  calculated  a  linear 
relationship  between  r  and  t500  at  each  station  from 
deep  hydrocasts  extending  to  the  bottom.  The  round 
trip  t  between  0  db  and  p  is 


where  p(p')  is  density,  g0  the  constant  acceleration  of 
gravity  at  9.8  m  s-2,  and  c(p')  the  sound  speed  (see 
Meinen  and  Watts  1997).  The  500-db  level  was  selected 
as  the  index  level  in  order  to  utilize  the  relatively  large 
number  of  hydrocasts  that  only  extended  to  this  depth. 
Because  the  waters  in  the  JES  below  500  db  are  nearly 
homogeneous,  a  simple  linear  relationship  together 
with  the  CTDs  shifts  the  measured  tp  records  that 
are  collected  at  various  depths  p  to  a  calibrated  t500 
[similar  to  the  methods  employed  by  Tracey  et  al. 
(1997)  and  Meinen  and  Watts  (1998)].  The  x500  accu- 
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Fig.  1.  The  Japan/East  Sea.  Solid  triangles  indicate  PIES  sites  and  white  circles  are  in  situ  CTD  surveys  by  KORDI  and  KODC. 

Bathymetry  contours  are  in  meters. 


racy  is  about  1.1  ms  when  compared  against  several 
CTDs.  Finally,  the  r500  records  were  OI  mapped  on  a 
0.125°  X  0.125°  grid  in  the  UB  using  a  42-km  correla¬ 
tion  length  scale,  as  was  used  by  Mitchell  et  al.  (2004). 

The  TRMM  Microwave  Imager  (TMI)  radiometer 
with  a  10.7-GHz  channel  collects  SST  data  through 
nonprecipitating  clouds,  which  is  a  distinct  advantage 
over  the  traditional  infrared  SST  observations  that  re¬ 
quire  a  cloud-free  field  of  view.  The  TRMM  SST  data 
were  downloaded  from  the  Remote  Sensing  Systems 
Web  site  as  daily  3-day  averages  and  weekly  composites 
on  a  0.25°  X  0.25°  grid  covering  the  UB.  After  filling  in 
missing  values  of  SST  by  interpolation,  we  produced 
736  daily  time  series  of  Ol-mapped  SST  on  the  same 
grid  as  Ol-mapped  tsoo  using  the  MOD  AS  climatologi¬ 
cal  SST  as  the  mean  field  in  the  OI  calculation. 


3.  MI-GEM  technique 

a.  Development  of  MI-GEM  fields 

Mitchell  et  al.  (2004)  illustrated  clearly  the  difficul¬ 
ties  in  constructing  a  GEM  field  in  the  UB.  They 
showed  two  pairs  of  temperature  profiles,  for  example 
(Fig.  2).  One  pair  was  collected  at  a  single  location  in 
the  summer-  and  wintertime,  and  the  other  pair  was 
collected  at  different  locations  in  the  same  month.  Each 
profile  pair  has  essentially  the  same  t500,  but  the  ele¬ 


ments  have  different  vertical  shapes,  for  example,  a 
warm  surface  layer  with  a  shallow  thermocline  versus  a 
mixed  layer  with  a  depressed  thermocline.  The  stan¬ 
dard  GEM  technique  estimates  a  single  T  profile  from 
t500  and,  hence,  the  standard  GEM  fields  cannot  distin¬ 
guish  between  these  profile  pairs.  We  utilize  available 
SST  data  to  improve  the  estimates.  Mitchell  et  al. 
(2004)  also  used  the  SST  data  in  their  residual  GEM 
technique,  but  they  classified  the  SST  into  only  three 
broad  bins.  In  contrast,  we  selected  SST  bins  having 
1°C  intervals  from  5°  to  30°C  in  order  to  index  with  a 
well-resolved  SST.  In  addition,  our  technique  estimates 
the  ML  structure  as  described  below.  As  Mitchell  et  al. 
(2004)  mentioned,  the  incorrect  estimation  of  the  ML 
was  the  largest  source  of  error  in  the  residual  GEM. 
Altogether,  our  MI-GEM  technique  depends  on  the 
presence  or  absence  of  an  ML,  and  utilizes  multi-index 
lookup  tables  of  T  and  8  as  functions  of  t500,  SST, 
and  p. 

Figure  3  illustrates  the  difference  between  the  stan¬ 
dard  GEM  and  MI-GEM  techniques,  using  scatterplots 
of  T  versus  t500.  The  standard  two-index  GEM  tech¬ 
nique  (e.g.,  Watts  et  al.  2001)  seeks  a  smooth,  single¬ 
valued  relationship  between  the  data  at  each  p  and  x500, 
exemplified  by  the  black  curve  on  each  scatterplot  at 
100,  200,  and  300  db.  The  superimposed  black  curves 
are  smoothing  cubic-spline  curves  that  are  fitted  to  all 
of  the  T  data  from  the  historical  hydrocasts,  as  repre- 
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Fig.  2.  (left)  Hydrographic  temperature  profiles  near  36.36°N,  134.08°E  for  Mar  and  Aug, 
illustrating  the  seasonally  changing  vertical  structure,  (right)  Hydrographic  temperature  pro¬ 
files  for  two  locations  that  are  325  km  apart  during  the  month  of  Feb.  In  each  panel,  two 
measured  temperature  profiles  (solid  lines)  with  similar  t500  values  are  compared  with  profiles 
estimated  from  the  standard  GEM  technique  (dashed  line)  (reproduced  from  Mitchell  et  al. 
2004). 


sented  by  black  dots.  The  standard  2D  GEM  lookup 
field  of  T,  as  a  function  of  tsoo  and  p,  is  produced  by 
fitting  a  cubic-spline  curve  to  the  T  versus  r500  data  on 
each  10-db  pressure  level  between  0  and  500  db,  gen¬ 
erating  51  spline  curves. 

The  MI-GEM  technique  regroups  the  same  T  pro¬ 
files  into  subsets  that  depend  on  SST  and  the  presence 
or  absence  of  a  mixed  layer.  As  an  illustration  for  just 
one  of  the  SST  bins  (15.5°-16.5°C),  the  red  and  blue 
curves  in  Fig.  3  are  fitted  to  the  subsets  of  NML  and  ML 
groups.  When  the  temperature  difference  between  0 
and  30  db  is  larger  (smaller)  than  1°C,  the  profile  was 
classified  as  being  an  NML  (ML)  group.  At  the  100-db 
level,  the  ML  group  has  higher  T  than  that  of  the  NML 
group  when  t500  is  less  than  0.676  s.  The  reverse  is  the 
case  for  the  200-db  level.  As  t500  decreases,  the  absolute 
differences  of  T  between  the  NML  and  ML  groups  in¬ 
crease  at  100  and  200  db.  At  the  300-db  level,  the  data- 
fitted  curves  for  both  of  the  groups  are  consistent  ev¬ 
erywhere  with  the  black  curve.  This  relationship  im¬ 
plies  that  a  warm  (cold)  temperature  in  the  UB  at  one 
layer  is  compensated  by  a  cold  (warm)  temperature  at 
the  other  layer  to  retain  a  selected  t500  value.  The  re¬ 


sults  of  this  regrouping  to  utilize  SST  and  ML  informa¬ 
tion  are  shown  in  Fig.  4.  The  profiles  were  first  grouped 
into  NML  and  ML  categories  and  sorted  according  to 
SST.  Then,  within  1°C  intervals  of  SST,  cubic-spline 
curves  were  fitted  on  each  pressure  level  at  10-db  in¬ 
tervals  to  generate  a  standard  GEM  field  of  T  for  each 
separate  group  as  a  function  of  r500  and  p.  Last,  it  was 
necessary  to  fit  smoothing  cubic  splines  along  the  SST 
axis  at  each  pressure  level  to  make  smoothly  varying 
MI-GEM  fields  (Fig.  4).  The  same  procedures  were 
applied  to  construct  MI-GEM  fields  of  8. 

To  interpret  PIES  data  using  this  new  technique,  we 
further  need  a  criterion  for  switching  between  the  NML 
and  ML  MI-GEM  fields.  The  histogram  in  Fig.  5  ex¬ 
hibits  the  monthly  distribution  of  the  number  of  NML 
and  ML  profiles  from  the  2248  hydrocasts  that  are  used 
to  construct  the  NML  and  ML  MI-GEMs.  This  histo¬ 
gram  illustrates  that  NML  and  ML  profiles  predomi¬ 
nantly  occur  during  the  summer  and  winter  months, 
respectively.  The  2-yr  time  series  of  NOGAPS  basin- 
averaged  wind  stress  through  the  UB  (Fig.  5a)  illus¬ 
trates  the  significant  seasonal  variation,  that  is,  strong 
(weak)  wind  during  the  winter  (summer)-time,  which 
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Fig.  3.  Scatterplots  of  T  vs  r500  at  (a)  100,  (b)  200,  and  (c)  300  db.  The  superimposed  black  lines  are 
the  smoothed  standard  GEM  curves  fitted  to  all  the  2248  historical  hydrocasts  (black  dots).  The  super¬ 
imposed  red  and  blue  lines  are  MI-GEM  curves  fitted  to  NML  and  ML  hydrocasts  (red  crosses  and  blue 
asterisks)  with  SST  between  15.5°  and  16.5°C. 


agrees  with  the  monthly  distribution  of  the  number  of 
NML  and  ML  profiles.  We  choose  10  May  as  the  date 
to  switch  from  the  ML  to  the  NML  MI-GEM  fields. 
Although  the  ratio  of  NML  to  ML  profiles  is  about 
50:50  in  September,  the  date  to  switch  back  from  the 
NML  to  ML  MI-GEM  fields  is  chosen  as  10  October, 
because  within  the  UB  most  of  the  ML  profiles  are  in 
this  month.  In  addition,  we  included  a  30-day  transi¬ 
tional  period  that  is  centered  on  each  of  the  two  switch¬ 
ing  dates,  during  which  time  weighting  functions  for  the 
two  MI-GEM  fields  linearly  increased  and  decreased 
(Fig.  5c),  in  order  to  produce  smoothly  varying  profiles 
around  the  switching  dates. 

The  Ol-mapped  SST  data  from  the  TRMM  satellite, 
which  are  the  fields  we  used  as  input  to  the  MI-GEM, 
are  compared  with  the  SST  data  from  the  in  situ  CTD 
casts  that  are  collected  by  KORDI  and  KODC.  The 
TRMM  SST  data  are  interpolated  to  coincide  with  each 
CTD  location.  Figure  6  exhibits  a  scatterplot  of  CTD 
SST  versus  TRMM  SST  and  a  least  squares  fitted  line 
that  is  superimposed  on  the  plot.  The  two  SST  mea¬ 
surements  are  highly  correlated  with  ±0.957(°C)  error 
bounds  (gray-shaded  zone)  that  contain  the  central 


50%  of  the  points.  Figure  6  demonstrates  that  the 
TRMM  SST  data  are  biased  slightly  high  when  the  SST 
is  lower  than  ~15°C.  Therefore,  we  corrected  the 
TRMM  SST  data  using  the  linear  relationship  as 
{TRMM  SST  -  1.291)/0.947. 

b.  Percent  variance  captured  by  MI-GEM  fields 

The  percent  variance  that  is  captured  by  the  MI- 
GEM  fields  for  T  and  8  as  a  function  of  pressure  is 
calculated  for  each  of  the  two  variables  by  1  -  afes,(p)l 
°?ot(p)>  where  ofes(p)  is  the  variance  of  residuals  of  the 
2248  CTDs  from  the  MI-GEM  fields,  and  afot  is  the 
total  variance  among  the  2248  CTDs  (Watts  et  al. 
2001).  Figure  7  compares  the  percent  variance  that  is 
captured  by  the  residual  GEM  and  MI-GEM  fields. 
Averaged  through  the  depth  range  of  0-200  db,  more 
than  92%  and  88%  of  the  variances  in  the  observed  T 
and  8  fields  are  captured  by  the  MI-GEM.  These  per¬ 
cent  variances  are  3%  and  4%  higher,  reducing  by 
about  one-third  the  unexplained  error  variance,  than 
those  of  the  residual  GEM.  Averaged  through  the 
depth  range  of  200-400  db,  MI-GEM  captures  73%  and 
67%  of  the  variances  for  the  two  variables,  which  is  an 
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Fig.  4.  MI-GEM  T  fields  for  (a)  non-mixed-layer  and  (b)  mixed-layer  groups.  The  color 

bar  is  temperature  (°C). 


improvement  over  the  residual  GEM,  which  captures 
68%  and  60%,  respectively. 

4.  Comparisons  with  observations 

The  441  T  profiles  from  the  in  situ  CTD  casts  (Fig.  1) 
that  are  collected  by  KORDI  (305  casts)  and  KODC 


(136  casts)  are  compared  against  profiles  that  are  esti¬ 
mated  from  both  the  residual  GEM  and  the  MI-GEM. 
This  is  an  independent  comparison  because  those  CTD 
casts  were  not  used  to  construct  the  residual  GEM  and 
MI-GEM  fields.  The  daily  OI  maps  of  t500  and  SST, 
described  in  section  2,  were  combined  to  estimate  T 
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(c)  Weighting  Function 


Fig.  5.  (a)  Magnitude  of  basin-averaged  wind  stress  in  the  UB.  The  90-day  low-pass-filtered  results  are 
superimposed  with  red  and  blue  curves,  representing  NML  and  ML  periods,  respectively,  (b)  Monthly 
distribution  of  NML  (red)  and  ML  (blue)  groups  from  historical  hydrocasts,  (c)  Weighting  function 
between  NML  (red)  and  ML  (blue)  groups  for  MI-GEM  application  during  the  2  yr  of  observations. 


profiles  on  a  0.125°  X  0.125°  grid  using  the  residual 
GEM  and  the  MI-GEM  lookup  tables  (Ol-mapped 
SST  errors  are  simulated  in  the  appendix).  These  two 
sets  of  T  profiles  were  then  interpolated  horizontally  to 
coincide  with  each  CTD  location.  The  difference  be¬ 
tween  the  estimated  and  actual  T  was  calculated  at  each 
10-db  level  for  all  441  CTD  casts,  and  the  rms  residual 
is  plotted  in  Fig.  8.  The  comparison  reveals  that  PIES- 
derived  T  estimates  using  the  MI-GEM  have  lower  rms 
differences  from  the  actual  T  profiles,  particularly  in 
the  30-90-  and  150-250-db-depth  ranges. 

The  statistical  improvement  in  Fig.  8  is  small,  because 
many  of  the  441  CTD  profiles  are  well  estimated  by 
either  technique.  However,  the  improvement  that  is  re¬ 
alized  by  the  MI-GEM  is  significant  at  20%-30%  of  the 
CTD  stations.  More  importantly,  it  represents  a  quali¬ 
tative  improvement  by  eliminating  false  T  inversions 


and  thermocline  structures,  as  illustrated  in  the  next 
two  figures. 

Figure  9  compares  vertical  T  sections  from  CTD  casts 
along  KODC  lines  106  and  104,  with  coinciding  (in 
space  and  time)  PIES-derived  T  sections  calculated  us¬ 
ing  the  residual  GEM  and  MI-GEM.  The  KODC  ob¬ 
servational  surveys  were  conducted  bimonthly  and, 
hence,  we  have  12  comparisons  of  the  T  sections  at  each 
KODC  line.  Figure  9  shows  six  examples  of  T  sections 
at  two  KODC  lines,  which  exhibit  the  improvements  of 
the  MI-GEM  estimates  over  the  residual  GEM. 

The  North  Korea  Cold  Current  and  the  East  Korea 
Warm  Current  usually  meet  and  form  the  subpolar 
front  around  38°-40°N.  KODC  line  106  (37.895°N)  is 
near  the  southern  boundary  of  the  subpolar  front  and 
has  been  reported  to  cross  the  path  of  the  ESIW  (Kim 
and  Kim  1983).  The  T  sections  illustrate  that  the  re- 
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Fig.  6.  Scatterplot  of  SST  from  in  situ  CTD  vs  SST  from  the 
TRMM  satellite.  The  slope  of  a  line  fitted  by  least  squares  is  0.947 
(thick  line)  and  the  gray-shaded  zone  indicates  error  bounds  that 
contain  at  least  50%  of  the  predictions. 

sidual  GEM  underestimates  the  depth  of  the  1°C  water 
surface  (Zx),  the  lower  boundary  of  the  ESIW,  at  line 
106,  except  during  February  2000  (Fig.  9,  left  panels).  In 
the  August  1999  residual  GEM  T  section,  the  water 
layer  between  1°  and  2°C  (Zl2)  is  unrealistically  thin 
around  129.5°E.  During  August  and  October  2000,  the 
ZX2  layer  is  too  narrow  along  the  whole  transect.  The 
April  and  June  2001  residual  GEM  T  sections  have  this 
same  problem  at  the  outer  boundaries  of  the  Ulleung 
warm  eddy  (UWE).  All  of  these  cases  of  narrow  Zl,2 
occurred  in  regions  where  the  thermocline  is  shallow. 
In  contrast,  the  MI-GEM  T  sections  exhibit  a  more 
realistic  Zx  2  layer  thickness  in  all  cases  and  agree  well 
with  the  CTD  observations  along  line  106. 

The  UWE  is  usually  centered  near  37°N,  131°E,  and 
line  104  spans  it  along  37°N.  In  the  summertime,  the 
residual  GEM  and  MI-GEM  T  sections  both  tend  to 
agree  with  observations,  as  exemplified  by  the  August 
1999  section.  In  the  wintertime,  the  surface  layer  of  the 
UWE  is  homogenized  by  strong  winds  and  the  mixed 
layer  deepens  beyond  100  db,  as  exemplified  in  four 
CTD  sections  from  December  1999,  March  and  April 
2000,  and  February  2001.  In  each  case,  the  residual 
GEM  T  sections  contain  unreasonable  T  inversions  be¬ 
tween  30  and  100  db.  The  residual  GEM  T  section  on 
line  106  for  February  2000  contained  the  same  fault  in 
the  upper  100  db.  In  contrast,  the  T  inversions  are 
eliminated  from  the  corresponding  MI-GEM  T  sections 
and  the  agreement  with  the  CTD-derived  ML  depth 
and  structure  is  qualitatively  better.  An  additional  im¬ 
provement  is  illustrated  in  the  June  2000  T  section,  in 
which  a  strongly  developed  lens-like  eddy  structure  that  • 


was  revealed  by  the  CTD  was  replicated  by  the  MI- 
GEM,  but  not  by  the  residual  GEM. 

Figure  10  compares  the  horizontal  maps  of  Z,  that 
are  calculated  from  CTD  surveys  and  PIES-derived 
maps,  interpreted  by  each  of  the  two  GEM  techniques 
for  the  same  months  as  the  KODC  line  106  sections 
shown  in  Fig.  9.  Along  the  northern  edge,  this  provides 
a  map  view  of  the  same  faults  with  the  residual  GEM 
estimates  that  are  “repaired”  by  the  MI-GEM  interpre¬ 
tation  of  the  PIES  data,  as  viewed  in  the  vertical  T 
sections  in  Fig.  9.  The  August  1999,  and  August  and 
October  2000  residual  GEM  Z,  maps  estimate  unreal¬ 
istically  shallow  Zx  near  38°N,  which  disagrees  with  the 
CTD  Zx  maps.  The  April  and  June  2001  residual  GEM 
maps  illustrate  excessively  steep  changes  of  Zx  at  the 
eastern  and  western  boundary  of  the  UWE.  The  MI- 
GEM  Zx  maps  eliminate  the  shallow  Zx,  and  agree  well 
with  the  CTD  Z,  maps  for  all  of  the  compared  cases. 

In  other  portions  of  these  maps  (Fig.  10)  two  other 
types  of  error  of  the  residual  GEM  estimates  are  found. 
In  the  February  2000  and  April  2001  maps,  the  spatial 
extent  of  the  UWE  in  the  residual  GEM  maps  is  too 
large  when  compared  with  that  in  the  CTD  and  MI- 
GEM  maps.  This  happens  because  the  residual  GEM 
overestimates  Zx  in  the  western  and  southern  parts  of 
the  UWE.  The  residual  GEM  estimates  also  err  signifi¬ 
cantly  in  the  southwestern  comer  of  the  mapping  do¬ 
main,  where  they  generate  a  feature  with  Zx  that  is 
deeper  than  400  m.  This  unrealistic  feature  does  not 
occur  in  the  CTD  and  MI-GEM  Zx  maps.  Figure  11 
demonstrates  that  this  erroneous  feature  originates 
from  a  bias  of  the  MODAS  dataset,  which  is  the  refer¬ 
ence  climatology  of  the  residual  GEM.  The  MODAS 
vertical  T  section  along  130°E  in  March  illustrates  this 
problem  south  of  36°N,  where  Zx  dives  unrealistically 
deep.  This  error  produced  the  isolated  deep  cold  water 
mass  that  is  estimated  by  the  residual  GEM. 

5.  Summary  and  discussion 

Our  objective  has  been  to  interpret  data  from  a  PIES 
array  that  was  deployed  for  2  yr  in  the  Ulleung  Basin  of 
the  southwestern  JES.  We  wanted  to  improve  upon 
estimates  of  T  and  8  profiles  that  were  obtained  using 
the  standard  GEM  technique,  because  variations  in  the 
shallow  permanent  thermocline  that  are  caused  by  ed¬ 
dies  are  difficult  to  separate  from  the  spatiotemporally 
varying  annual  signal.  Mitchell  et  al.  (2004)  addressed 
this  difficulty  by  developing  a  residual  GEM  technique. 
That  technique  demonstrated  considerable  success  in 
mapping  the  100-m  T  field,  which  has  served  as  an  ef¬ 
fective  proxy  for  mapping  currents  and  eddies  in  the 
JES.  However,  the  residual  GEM  technique  still  pro- 
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Fig.  7.  Percent  of  (a)  T  and  (b)  8  signals  captured  by  the  residual  GEM  (thin  line)  and  MI-GEM  (thick  line). 


duced  several  erroneous  features,  with  the  most  prob¬ 
lematic  being  significant  T  inversions  within  the  near¬ 
surface  layer  in  the  wintertime.  That  technique  also 
estimated  an  unreasonably  thin  intermediate  water 
layer  in  regions  where  the  thermocline  was  very  shal¬ 
low.  Unrealistically  thick  intermediate  water  layers 
were  also  obtained  in  one  region  where  the  MOD  AS 
climatology  profiles  have  bias  errors  as  demonstrated 
by  comparisons  with  concurrent  in  situ  CTD  observa¬ 
tions. 

From  our  subsequent  studies,  we  now  understand 
that  because  the  residual  GEM  technique  profiles  of 
the  T  and  S  anomaly  (relative  to  MOD  AS  climatology) 
have  subsurface  extrema,  when  these  are  associated 
with  a  different  residual  t  (t  anomaly)  and  added  back 


Temperature  (°C) 


Fig.  8.  The  rms  differences  of  T  between  441  CTD  profiles  and 
PIES-derived  profiles  estimated  by  the  residual  GEM  (thin  line) 
and  the  MI-GEM  (thick  line).  The  dashed-dotted  line  simulates 
added  SST  error,  as  explained  in  the  appendix. 


to  the  smooth  MODAS  climatological  profiles,  the 
summed  estimate  can  generate  false  inversions.  The 
subsurface  positive  (negative)  anomaly  extremum  that 
is  estimated  when  a  large  negative  (positive)  residual  t 
occurs  is  typically  associated  with  a  strong,  warm  (cold) 
eddy.  For  example,  during  the  wintertime,  the  UWE 
has  a  large  negative  t  anomaly.  The  residual  GEM  cor¬ 
rects  the  climatological  T  profiles  by  adding  a  positive 
subsurface  extremum,  which  cause  an  unreal  subsur¬ 
face  T  inversion,  as  examplified  in  Fig.  9.  Yet  none  of 
the  original  data  profiles  have  subsurface  T  inversions. 
So,  we  returned  to  seek  a  technique  that  developed 
GEM  fields  directly  from  the  original  data  profiles 
rather  than  anomaly  profiles.  However,  we  also  simul¬ 
taneously  addressed  the  shallow  pycnocline  problem 
(of  separating  the  eddy  signal  from  the  annual  signal) 
by  including  additional  information  that  was  available 
from  SST  and  the  temporal  presence  or  absence  of 
near-surface  mixed  layers  in  the  historical  hydrocasts. 
We  selected  the  SST  bins  having  1°C  intervals  from  5° 
to  30°C,  while  Mitchell  et  al.  (2004)  classified  the  re¬ 
sidual  SST  into  only  three  bins. 

We  call  this  new  technique  an  MI-GEM,  parameter¬ 
ized  by  t50o,  SST,  ML  information,  and  pressure  to 
make  multi-index  lookup  tables  for  T  and  8  from  his¬ 
torical  hydrocasts.  To  parameterize  the  ML  strength, 
the  historical  hydrocasts  are  regrouped  into  NML  and 
ML  groups.  The  appropriate  dates  for  NML  and  ML 
fields  are  determined  by  the  monthly  distribution  of  the 
number  of  NML  and  ML  profiles  in  the  historical  hy¬ 
drocasts,  which  also  correlated  with  the  strength  of  the 
90-day  low-pass-filtered  winds  during  our  observation. 
The  MI-GEM  technique  reduced  by  about  one-third 
the  unexplained  error  variance  of  the  residual  GEM. 
Although  the  MI-GEM  technique  includes  the  ML  in¬ 
formation,  it  is  still  only  an  estimate  of  the  ML  struc- 
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Fig.  9.  Vertical  T  sections  along  lines  (left)  106  and  (right)  104  from  KODC  bimonthly  surveys  of  serial  oceanographic 
stations  (red  dots)  compared  to  those  estimated  by  the  residual  GEM  and  the  MI-GEM  techniques.  Magenta  lines  in  line 
106  sections  indicate  the  Z12  layer  boundaries  (1°,  2°C).  The  color  bar  is  temperature  (°C),  and  the  contour  interval  is  1°C. 


ture,  and  this  remains  the  largest  source  of  unexplained 
error  variance. 

The  MI-GEM  technique  is  applied  to  the  t  measured 
by  a  PIES  array  in  the  UB  to  produce  water  column 


estimates  of  T  and  8  from  0  to  500  db.  At  first  sight,  the 
comparison  of  441  T  profiles  from  in  situ  CTD  casts 
against  coinciding  residual  GEM  and  MI-GEM  T  pro¬ 
files  reveals  only  small  a  quantitative  improvement. 
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KODC  stations 


Fig.  10.  Horizontal  maps  of  Z x  calculated  from  KODC  bi¬ 
monthly  surveys  (red  dots)  compared  to  those  from  the  residual 
GEM  and  the  MI-GEM  estimates.  The  color  bar  is  depth  (m),  and 
the  contour  interval  is  20  m. 


MODAS  (March,  along  130°E) 


FIg.  11.  A  vertical  T  section  of  MODAS  along  130°E  in  Mar, 
illustrating  a  problem  in  the  MODAS  T  field  south  of  36°N.  Ma¬ 
genta  lines  indicate  the  Zh2  layer  boundaries  (1°,  2°C). 


Moreover,  the  T  maps  at  a  100-m  depth  from  the  MI- 
GEM  estimates  are  almost  identical  to  those  from  the 
residual  GEM  estimates  (not  shown  in  this  paper). 
However,  the  MI-GEM  technique  eliminates  unrealis¬ 
tic  T  inversions  and  thin  or  thick  intermediate  water 
layers,  which  were  misestimated  by  the  residual  GEM. 
These  significant  qualitative  improvements  are  demon¬ 
strated  by  comparisons  of  vertical  T  sections  and  hori¬ 
zontal  Zj  maps  that  are  calculated  from  independent 
KODC  bimonthly  in  situ  CTD  surveys.  The  small  re¬ 
duction  in  scatter  around  the  MI-GEM  fields,  as  shown 
in  Fig.  5,  results  in  smaller  error  bars  on  the  estimated 
vertical  T  and  8  profiles  from  r  and,  hence,  the  MI- 
GEM  technique  almost  always  produced  improved  up¬ 
per  water  column  profiles  of  T  and  8  compared  to  the 
residual  GEM  technique.  The  qualitative  improvement 
of  the  intermediate  water  layer  estimates  provides  us  an 
opportunity  to  investigate  the  ESIW  dynamics  using 
daily  maps  of  the  ESIW  circulation  pattern  in  future 
work. 

For  the  above  reasons  the  MI-GEM  is  clearly  pre¬ 
ferred  over  the  residual  GEM  for  Ulleung  Basin  stud¬ 
ies.  Three  points  address  the  naturally  arising  question 
as  to  whether  the  MI-GEM  technique  should  outper¬ 
form  other  techniques  in  all  ocean  regions.  1)  The  ob¬ 
vious  point  is  that  multiple-index  proxies  should  always 
reduce  the  unexplained  variance  compared  to  the  stan¬ 
dard  two-index  GEM,  provided  that  one  selects  addi¬ 
tional  indices  that  independently  account  for  a  portion 
of  the  observed  variability.  2)  It  remains  possible  that  a 
residual  GEM  may  provide  better  estimates  of  vertical 
profiles  in  a  location  where  a  good  climatology  is  avail- 
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able  and  the  climatological  annual  signal  is  large  com¬ 
pared  to  the  mesoscale  eddy  signal.  Residual  GEM 
fields  might  also  be  improved  by  grouping  stations  as 
we  did  here.  Nevertheless,  the  residual  GEM  technique 
should  always  be  approached  cautiously,  because  the 
anomaly  fields  that  it  generates  may  contain  subsurface 
extrema,  which  are  intrinsically  subject  to  the  problem 
of  estimating  unrealistic  subsurface  temperature  and 
density  features.  3)  The  lookup  tables  that  are  gener¬ 
ated  by  the  MI-GEM  technique  (and  the  standard 
GEM)  are  composed  directly  from  real  hydrographic 
profiles.  Every  input  profile  is  stably  stratified  and  con¬ 
tributes  a  point  to  every  estimated  output  depth. 
Hence,  the  MI-GEM  smoothed  lookup  tables  should 
be  relatively  immune  to  containing  unreal  features. 
Moreover,  the  MI-GEM  fields  may  be  examined  di¬ 
rectly,  rather  than  as  anomaly  fields;  and  if  some  unre¬ 
alistic  feature  were  generated  it  would  come  to  a  user’s 
attention  readily. 

In  regions  where  the  seasonal  variations  are  suffi¬ 
ciently  decoupled  vertically  from  the  large-scale  pycno- 
cline  variations,  the  standard  GEM  (Meinen  and  Watts 

2000)  and/or  the  improved  standard  GEM  (Watts  et  al. 

2001)  will  probably  satisfactorily  capture  a  high  per¬ 
centage  of  the  variance  of  the  temperature  and  specific 
volume  anomaly  profiles.  In  all  regions  of  the  ocean,  it 
seems  likely  that  the  GEM  lookup  tables  can  be  im¬ 
proved  by  using  additional  proxy  indices  to  reduce  the 
residual  variance.  This  paper  has  demonstrated  the  par¬ 
ticular  advantages  of  combining  t  data  with  satellite 
SST  and  ML  information  in  an  MI-GEM. 
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APPENDIX 

Simulated  Introduction  of  SST  Errors 

To  evaluate  the  impact  of  errors  in  the  Ol-mapped 
TRMM  SST  data,  we  conducted  a  simulation  of  esti¬ 
mating  T  profiles  using  the  MI-GEM  technique  after 
adding  a  1°C  error  to  all  of  the  SST  data.  The  estimated 


T  profiles  were  interpolated  to  coincide  with  each  CTD 
location,  and  then  the  rms  residual  between  estimated 
and  actual  T  was  included  in  Fig.  8  (dashed-dotted 
line).  This  simulated  additional  SST  error  increased  the 
rms  residual  in  the  upper  100  m.  However,  even  if  the 
SST  data  had  a  1°C  error  added,  the  T  profiles  that 
were  estimated  from  the  MI-GEM  were  still  better  than 
those  estimated  from  residual  GEM  without  the  SST 
error. 
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